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ABSTRACT 


A revised  version  of  the  SHAKE  program  was  prepared  and  used  to  study  the  effects  of 
subsurface  conditions  on  the  earthquake  ground  motion  in  the  Loma  Prieta  earthquake. 
Preliminary  soil  profile  data  from  the  sites  of  the  Oakland  Outer  Harbor  Wharf  and  Apeel  2 
strong  motion  stations  are  used  to  calculate  ground  motions,  which  are  then  compared  with 
the  recorded  ground  motions  using  response  spectra  calculated  for  a 5%  damping  ratio. 
Parameters  affecting  the  amplitude  of  the  calculated  ground  motion  are  examined.  Response 
Spectra  for  recorded  and  calculated  ground  motions  are  compared  with  recommended  design 
spectra  (NEHRP,  1 988).  It  is  shown  that  for  periods  less  than  1 .4  s the  response  spectra  for 
recorded  far  source  earthquake  motions  at  Oakland  Wharf  and  Apeel  2 fall  outside  the 
envelope  of  the  applicable  design  spectra,  and  that  response  spectra  for  deeper  soil  profiles 
calculated  for  near  source  conditions  exceed  the  design  spectra  by  a considerable  margin. 
Design  spectra  for  the  San  Francisco  Bay  region,  recently  proposed  in  a USGS  study,  are 
reasonably  close  to  the  calculated  near  source  spectra  for  deeper  soil  profiles  for  periods  less 
than  1 s,  but  they  are  conservative  for  the  bedrock  motion,  and  extremely  conservative  for 
longer  period  structures. 

Key  Words:  dynamic  soil  properties;  earthquake  engineering;  ground  motion;  shear  wave 
propagation;  response  spectra;  soil  dynamics;  wave  mechanics 
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1. 


INTRODUCTION 


One  of  the  prominent  features  of  the  October  17,  1989  Loma  Prieta  Earthquake  was  the  close 
correlation  between  the  subsurface  conditions  and  the  damage  pattern.  While  this  is  not  a 
new  observation  (re:  Mexico  City,  Caracas,  etc.),  it  appears  in  this  instance  that  even  very 
stiff  structures  were  damaged  on  sites  with  relatively  deep  soil  deposits,  while  similar 
structures  in  the  vicinity,  resting  on  shallow,  competent  soil  or  rock,  suffered  no  damage. 
Substantial  amplification  of  the  earthquake  motion  was  observed  in  areas  of  Bay  Mud  deposits 
(average  amplification  of  maximum  vertical  and  horizontal  accelerations  by  factors  of  1 .8  and 
2.6,  respectively)  and  on  sites  of  deep  alluvial  deposits  (average  amplification  of  maximum 
vertical  and  horizontal  accelerations  by  factors  of  1.9  and  1.8,  respectively)  (EERI,  1989). 
Coupled  with  the  longer  dominant  period  of  the  horizontal  ground  motion  of  most  deep  soil 
deposits,  the  horizontal  acceleration  amplifications  result  in  much  larger  amplifications  of 
velocities  and  displacements. 

After  the  Loma  Prieta  earthquake,  the  National  Institute  of  Standards  and  Technology  (NIST) 
initiated  a research  effort  which  was  designed  to  address  gaps  in  the  state  of  the  art  of 
earthquake  engineering  and  deficiencies  in  design  standards  which  became  apparent  as  a 
result  of  the  earthquake  damage.  One  of  the  topics  identified  as  requiring  further  study  is  the 
effect  of  subsurface  conditions  on  the  earthquake  ground  motion.  The  research  undertaken 
by  NIST  focussed  on  the  prediction  of  the  effect  of  the  intervening  unconsolidated  deposits 
on  the  propagation  of  earthquake  ground  motions  from  the  underlying  bedrock  to  the  ground 
surface,  as  well  as  on  the  selection  of  bedrock  motions  for  seismic  design. 

Strong  motion  data  collected  during  the  earthquake,  coupled  with  information  acquired  in 
subsequent  subsurface  investigations,  afford  the  opportunity  to  test  the  validity  of  present 
mathematical  models,  dealing  with  the  propagation  of  the  earthquake  motion  from  the 
underlying  bedrock  to  the  ground  surface.  In  one  phase  of  the  NIST  research,  which  is  carried 
out  at  the  University  of  California,  Davis,  subsurface  data  from  the  site  of  the  strong  motion 
recording  on  Treasure  Island  are  utilized  to  examine  mathematical  models  for  predicting  the 
ground  motion.  In  this  latter  project  strong  motion  records  from  rock  outcrop  motions  at 
Yorba  Buena  and  Rincon  Hill  were  utilized  to  generate  a ground  surface  motion,  which  then 
could  be  compared  with  the  strong  motion  record  from  Treasure  Island. 

In  this  report,  a revised  version  of  the  SHAKE  program  (Schnabel  et  al.,  1972),  prepared  as 
part  of  the  NIST  research,  is  used  to  study  the  effect  of  various  soil  conditions  encountered 
in  the  San  Francisco  Bay  area  on  the  horizontal  components  of  the  ground  motion.  Vertical 
ground  motions  are  not  studied  in  this  NIST  project.  A user  manual  for  the  revised  SHAKE 
program  will  be  published  at  a later  date.  Some  of  the  revisions  are  discussed  in  the  following 
section. 
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2. 


MODELING  THE  PROPAGATION  OF  THE  BEDROCK  MOTION  TO  THE  GROUND 
SURFACE 


Plots  of  dynamic  stress-strain  properties,  used  in  this  study  for  Bay  Mud,  alluvium,  and  sand 
under  small  confining  pressures  are  shown  in  Figure  2.1 . The  shear  modulus  is  plotted  as  a 
function  of  G nulx,  the  shear  modulus  at  very  small  strains,  such  as  those  associated  with  shear- 
wave  propagation  velocity  measurements  in  a geophysical  soil  explorations.  Damping  is  given 
as  a percentage  of  critical  damping.  It  can  be  seen  that  the  dynamic  shear  modulus 
decreases,  and  the  damping  ratio  increases  with  increasing  strains. 

Many  different  methods  can  be  used  to  model  the  propagation  of  the  bedrock  motion  to  the 
ground  surface,  such  as  modeling  of  wave  propagation,  lumped  mass  models,  and  finite 
element  models.  An  important  aspect  of  a numerical  model  is  whether  or  not  it  can 
reasonably  represent  the  non  linear  behavior  of  soils  under  cyclic  loading  conditions.  From 
this  viewpoint  the  models  fall  into  two  categories:  equivalent  linear,  and  non  linear  models. 
In  the  equivalent  linear  models,  the  dynamic  shear  modulus  and  the  damping  ratio  for  the 
entire  time  history  of  the  motion  are  fixed  for  each  soil  layer  for  a specific  strain  level  which 
is  defined  as  a fraction  of  the  maximum  strain  experienced.  This  approach  does  not  permit 
consideration  of  cyclic  degradation  effects  on  stiffness  and  damping.  It  is  therefore  not 
sensitive  to  the  number  of  strain  cycles  applied,  although  allowance  for  the  magnitude  of  the 
earthquake  can  be  made  by  the  setting  of  the  user  defined  strain  level  for  which  the 
equivalent  linear  properties  are  determined.  The  equivalent  linear  approach  can  also  not  be 
used  to  predict  the  ground  motion  for  sites  which  experience  cyclic  mobility  or  liquefaction. 
However,  it  is  extensively  used  to  predict  the  liquefaction  potential  of  sites  in  terms  of  the 
cyclic  stress  or  strain  level  before  liquefaction  occurs.  Another  limitation  to  the  use  of 
equivalent  linear  models  is  the  magnitude  of  cyclic  strains  experienced  in  the  earthquake.  If 
these  strains  are  very  large,  most  soils  will  experience  cyclic  degradation  and  the  use  of  an 
equivalent  linear  model  is  therefore  inappropriate. 

In  this  study  the  SHAKE  program  (Schnabel  et  al.,  1972)  was  used,  which  is  an  equivalent 
linear  program,  based  on  one-dimensional  wave  propagation  theory.  The  1-D  wave  equation 
is  used  on  the  basis  of  the  assumption  that  the  soil  deposit  is  horizontally  layered 
(orthotropic).  The  ground  motion  is  idealized  as  in-plane  horizontal  displacements,  propagated 
upward  from  an  isotropic  halfspace  of  bedrock  with  a horizontal  surface.  This  idealization 
is  reasonably  applicable  to  many  situations,  particularly  in  sedimentary  deposits.  Equations 
used  in  the  SHAKE  program,  which  were  derived  by  Schnabel  et  al.,  1972,  and  amended  by 
Udaka  and  Lysmer,  1973,  are  reviewed  in  the  appendix.  Since  it  was  authored  in  1972,  the 
SHAKE  program  has  undergone  successive  revisions.  The  last  version  published  before 
initiation  of  the  NIST  ground  motion  study  was  by  Sun  and  Galesorkhi,  1988.  In  this  latter 
version,  the  number  of  soil  types  for  which  dynamic  properties  can  be  specified  was  increased 
from  4 to  13,  and  specific  dynamic  properties  for  clays  and  sands  were  proposed.  In  the 
current  revision  the  number  of  soil  layers  in  the  profile  for  which  the  calculations  can  be 
performed  was  increased  from  20  to  50.  This  not  only  enables  us  to  better  represent  complex 
subsurface  conditions,  but  also  permits  the  use  of  thinner  layers.  It  was  reasoned,  that  in 
order  to  properly  represent  the  higher  frequency  range  of  the  ground  motion,  layer  thicknesses 
should  not  exceed  1/4  of  the  wave  length.  Other  revisions  in  the  program  are  discussed  in 
more  detail  in  the  user  manual  for  the  revised  program  (Idriss  et  al.,1992). 
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Figure  2.1:  Dynamic  Stress-Strain  Properties  of  Soils  Used  in  this  Study 
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3.  SOIL  PROFILES  STUDIED 

3.1  Soil  Properties  Used 

3.1.1  Choice  of 

Figure  1 shows  the  dynamic  shear  modulus  reduction  G/G ^ as  a function  of  strain,  relative 
to  the  dynamic  shear  modulus  at  very  small  strains,  G^.  This  still  leaves  the  problem  of 
determining  a value  for  G^.  In  practice  this  is  usually  accomplished  in  situ  by  measuring  the 
shear  wave  velocity.  The  shear  wave  velocity  is  related  to  G^  by  the  following  relationship: 


. . . (1) 


max 


where : p = mass  density  of  soil 

Vs  = shear  wave  propagation  velocity 


In  sands,  the  shear  modulus  has  been  found  to  be  approximately  proportional  to  the  square 
root  of  the  mean  normal  effective  stress  (Seed  & Idriss,  1970): 


. . . (2) 


where:  o0  = mean  normal  effective  stress 


K2  — empirical  constant 
G = shear  modulus 


These  values  can  be  converted  to  SI  units  as  follows: 


• . . (3) 


where:  o0  = mean  normal  stress  in  MPa 

K^OcPa)05  = .2188-K2(psf)05  = empirical  constant 
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The  following  soil  properties  were  used  in  the  analysis  of  the  San  Francisco  area  profiles. 
These  properties  represent  an  estimate  of  average  values.  For  example,  shear  wave  velocity 
values  as  low  as  55  m/s  (180  ft/s)  have  been  reported  for  some  New  Bay  Mud  deposits. 

For  fill: 

Kxmat)  = 50  (psf)as;  or  = 10.94  (kPa)05 

Unit  Weight:  .1  Kef  —1.6  Mg/m3 

For  Alluvium: 

K2(max)  = 60  to  90  (psff5;  or  = 13  to  20  (kPaf5 

Unit  Weight  = .13  to  .14  Kef  = 2 to  2.2  Mg/m3 

For  New  Bay  Mud: 

Vs  = 300  to  600 ft/s  = 90  to  180  m/s 

Unit  Weight:  .11  to  .13  Kef  = 1.76  to  2.1  Mg/m3 

For  Old  Bay  Mud: 

Vs  = 1200  to  2500 ft/s  (for  600 ft  deep  deposit)  = 365  to  762  m/s 
Unit  weight  = .13  Kef  = 2.08  Mg/m3 

For  the  study  of  the  generic  soil  profiles,  representing  generic  soil  profile  types  specified  in  the 
National  Earthquake  Hazard  Reduction  Program  (NEHRP)  recommended  provisions  (NEHRP, 
1 988),  shear  wave  velocities  of  soft  soil  deposits  were  assumed  to  vary  from  200  to  400  ft/s 
(60  to  120  m/s). 

3.1.2  Relationship  between  and  fi  and  cyclic  shear  strain 

The  variations  of  modulus  reduction  ( G/G and  damping  ratio  (fi)  with  cyclic  strains  used 
in  this  study  are  shown  in  figure  2.1 . The  following  data  were  used  in  the  analysis: 
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Curves  NO  1 and  2 were  used  for  modulus  and  damping,  respectively,  of  New  and  Old  Bay 
Mud: 

CURVE  NO.  1:  modulus  for  clay  (Seed  & Sun,  1989),  upper  range. 

CURVE  NO.  2:  damping  for  clay  (Idriss,  1990). 


CURVE  NO. 

1 

CURVE  NO. 

2 

STRAIN 

G/Gmax 

STRAIN 

DAMPING 

.1 00E-03 

1.000 

.1 00E-03 

.24 

.300E-03 

1.00 

.300E-03 

.42 

.1 00E-02 

1.000 

.100E-02 

.80 

.300E-02 

.981 

.300E-02 

1.40 

.100E-01 

.941 

.1 00E-01 

2.80 

.300E-01 

.847 

.300E-01 

5.10 

.100E  + 00 

.656 

.100E  + 00 

9.80 

.300E  + 00 

.438 

.300E  + 00 

15.50 

.100E  + 01 

.238 

.100E  + 01 

21.00 

.300E  + 01 

.170 

.300E  + 01 

25.00 

.100E  + 02 

.130 

.100E  + 02 

28.00 

Curves  NO  3 and  4 were  used  for  modulus  and  damping,  respectively,  of  Alluvium. 


CURVE  NO.  3:  modulus  for  sand  (Seed  & Idriss,  1970)  - upper  Range. 

CURVE  NO.  4:  damping  for  sand  (Idriss  1990)  - (approximately  Lower  Range  from  Seed  and 
Idriss,  1970,  and  identical  to  damping  curve  used  for  Bay  Mud). 


CURVE  NO. 

3 

CURVE  NO. 

4 

STRAIN 

G/G^ 

STRAIN 

DAMPING 

.1 00E-03 

1.000 

.1 00E-03 

.24 

.300E-03 

1.000 

.300E-03 

.42 

.1 00E-02 

.990 

.1 00E-02 

.80 

.300E-02 

.960 

.300E-02 

1.40 

.1 00E-01 

.850 

.1 00E-01 

2.80 

.300E-01 

.640 

.300E-01 

5.10 

.100E  + 00 

.370 

.100E  + 00 

9.80 

.300E  + 00 

.180 

.300E  + 00 

15.50 

.100E  + 01 

.080 

.100E  + 01 

21.00 

.300E  + 01 

.050 

.300E  + 01 

25.00 

.100E  + 02 

.035 

.100E  + 02 

28.00 
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Curves  NO  17  and  18  were  used  for  modulus  and  damping,  respectively,  of  fill. 

CURVE  NO.  17:  sand,  average  confining  pressure  < 100  kPa  (Sun,  1988). 
CURVE  NO.  18:  damping  of  sand,  (Seed  & Idriss,  1970). 


CURVE  NO. 17  CURVE  NO. 18 


STRAIN 

G/Gmax 

STRAIN 

DAMPING 

.100E-03 

1.000 

.1 00E-03 

1.00 

.31  6E-03 

.978 

.100E-02 

1.60 

.100E-02 

.934 

.300E-02 

3.12 

.31 6E-02 

.838 

.1 00E-01 

5.80 

.100E-01 

.672 

.300E-01 

9.50 

.31 6E-01 

.463 

.100E  + 00 

15.40 

.100E  + 0Q 

.253 

.300E  + 00 

20.90 

.31 6E  + 00 

.140 

.100E  + 01 

25.00 

. 1 0QE  + 01 

.057 

.100E  + 02 

25.50 

3.2  Soil  Profiles  Studied 

The  soil  profiles  used  in  this  study  are  outlined  in  Table  1.  Profiles  No.  14  and  15  were 
obtained  from  preliminary  field  data  from  the  sites  of  the  Oakland  Outer  Harbor  Wharf  and 
Apeel  2,  a United  States  Geological  Service  (USGS)  array  in  Redwood  City,  CA,  where  a 
ground  motion  time  history  is  available.  These  latter  profiles,  as  interpreted  in  this  study  for 
input  into  the  SHAKE  program,  are  shown  in  figures  3.1  and  3.2. 
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TABLE  1:  SOIL  PROFILES 


Layer  Thickness,  m 


Depth  to 

New  Bay  Mud  Alluvium  Old  Bay  Mud  Water  Table 


3 

6 

15 

30 

0 

0 

0 

3 

30 

30 

3 

3 

6 


3 

3 

3 

3 

3 

12 

24 

24 

24 

3 

30 

3 

6 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

30 

150 


OAKLAND  WHARF  (figure  2) 
A P E E L 2 (figure  3) 

9 (V ’ = 60  m/s) 


9 (Vs  = 120  m/s) 
21  (Vs  = 60  m/s) 
21  (Vs  = 120  m/s) 
46  (Vy  = 60  m/s) 
46  (Vs  = 120  m/s) 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 
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OAKLAND  OUTER  HARBOR  WHARF.  Preliminary  Soil  Data 


Layer 

Soil  Type 

Depth 

Thickns. 

Unit  Wt. 

v. 

ft 

ft 

kef 

ft/s 

1 

Sand 

0 

6 

0.103 

440 

2 

Fine  Sand 

6 

6 

0.105 

440 

3 

Fine  Sand 

12 

6 

0.105 

440 

4 

Gray  Clay 

18 

4 

0.110 

440 

5 

Fine  Sand 

22 

5 

0.110 

440 

6 

0 m 

27 

10 

0.110 

920 

7 

0 0 

37 

10 

0.110 

920 

8 

0 

47 

10 

0.120 

1,280 

9 

Fine  Sand 

57 

9 

0.120 

1,280 

10 

Gray  Clay 

66 

6 

0.111 

770 

11 

Stiff  Clay 

72 

4 

0.120 

770 

12 

Gray  Clay 

76 

6 

0.113 

770 

13 

Sand 

82 

5 

0.120 

770 

14 

Clay 

87 

6 

0.120 

770 

15 

• 

93 

10 

0. 120 

770 

16 

" 

103 

10 

0.120 

770 

17 

• 

113 

10 

0.120 

770 

18 

m 

123 

10 

0.130 

1,070 

19 

Clay 

133 

10 

0.130 

1,070 

20 

Sand 

143 

6 

0. 120 

1,070 

21 

Silty  Clay 

149 

7 

0.130 

1,070 

22 

Silty  Clay 

156 

7 

0.130 

1,070 

23 

Gray  Sand 

163 

3 

0.125 

1,070 

24 

Clay 

166 

20 

0.130 

1,070 

25 

Clay 

186 

21 

0.130 

1,070 

26 

Sand 

207 

3 

0.120 

1,070 

27 

Clay 

210 

6 

0.130 

780 

28 

" 

216 

10 

0. 130 

780 

29 

Clay 

226 

10 

0. 130 

780 

30 

Sand 

236 

2 

0.125 

780 

31 

Clay 

238 

12 

0.130 

2,060 

32 

Sand 

250 

4 

0.125 

2,060 

33 

Clay 

254 

6 

0.130 

2,060 

34 

■ 

260 

10 

0.130 

780 

35 

* 

270 

10 

0.130 

780 

36 

• 

280 

10 

0.130 

780 

37 

• 

290 

12 

0.130 

1,400 

38 

Sand 

302 

22 

0.125 

1,400 

39 

Stiff  Clay 

324 

6 

0.130 

1,400 

40 

Stiff  Clay 

330 

10 

0.130 

1,400 

41 

Sandy  Loam 

340 

10 

0.125 

1,400 

42 

Sandy  Loam 

350 

11 

0.130 

1,400 

43 

Clay 

361 

11 

0.130 

1,400 

44 

Sand 

372 

12 

0.125 

1,400 

45 

Clay 

384 

10 

0.130 

1,400 

46 

Clay 

394 

10 

0.130 

1,400 

47 

Sand 

404 

20 

0. 130 

1,400 

48 

Sand 

424 

20 

0. 130 

1,400 

49 

Sand 

444 

50 

0.130 

2,380 

50 

Bedrock 

494 

Vs  = measured  shear  wave  velocity;  Unit  Weight  assumed; 

1 ft  = 0.3048  m;  1 kef  = 16.03  Mg/m3 

Figure  3.1:  Soil  Profile  near  California  Department  of  Mines  and  Geology  (CDMG)  Strong 
Motion  Station  No.  58472  at  Oakland  Outer  Harbor  Wharf. 
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APEEL  2 (Redwood  City),  Preliminary  Soil  Data 


Layer 

Soil  Type 

Depth  Thickns. 

Unit  Wt. 

Vs 

ft 

ft 

kef 

ft/s 

1 

Sdy.Cl.Fill 

3.61 

3.61 

.103 

574.4 

2 

Silty  Clay 

2.95 

.110 

180.4 

3 

(Lt.  Gray, 

3.28 

.1 1 1 

180.4 

4 

soft) 

3.28 

.1 1 1 

180.4 

5 

it 

3.28 

.112 

180.4 

6 

i» 

3.28 

.113 

180.4 

7 

n 

3.28 

.114 

180.4 

8 

ii 

3.28 

.114 

180.4 

9 

n 

31.2 

4.92 

.115 

180.4 

10 

Clay, (Yellowish 

4.92 

.13 

639.76 

1 1 

Brown,  Hard) 

6.56 

.13 

639.76 

12 

ii 

6.56 

.13 

639.76 

13 

ii 

6.56 

.13 

639.76 

14 

n 

6.56 

.13 

639.76 

15 

if 

65.9 

3.61 

.13 

639.76 

16 

Silty  Clay, 

6.23 

.13 

574.15 

17 

(Gray  w.  Olive 

6.23 

.13 

574.15 

18 

Mottles,  Stiff) 

6.23 

.13 

574.15 

19 

•i 

6.23 

.13 

574.15 

20 

it 

100.7 

9.84 

.13 

574.15 

21  Dense  Deposits  (no  data) 

Vs  = measured  shear  wave  velocity 
Unit  Weight  was  assumed 
1 ft  = 0.3048  m;  1 kef  = 16.03  Mg/m3 


Figure  3.2:  Soil  Profile  near  Strong  Motion  Station  USGS  Apeel  2,  Redwood  City,  CA. 
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4. 


RESULT  OF  RESPONSE  CALCULATIONS 


4.1  Presentation  of  results 

Since  the  calculations  produce  a great  number  of  data  points,  the  most  practical  way  to  study 
the  results  is  from  graphical  presentations.  What  we  are  comparing  is  the  rock  outcrop 
motion  used  to  calculate  the  bedrock  motion  at  the  base  of  the  deposit  and  the  calculated 
ground  motion  at  the  surface  of  the  deposit.  Since  ground  motions  are  measured  in  the  field 
as  acceleration  time  histories,  they  are  presented  in  this  form  in  the  report.  However, 
acceleration  time  histories,  even  when  plotted,  are  difficult  to  compare.  One  parameter  that 
can  be  compared  is  the  maximum  acceleration  in  the  record.  However,  such  a comparison 
is  insufficient,  because  this  quantity  alone  does  not  characterize  the  entire  acceleration  time 
history. 

One  effective  way  of  comparing  ground  motions  is  by  comparing  their  effect  on  structures 
with  different  fundamental  periods.  This  is  done  by  means  of  response  spectra,  which  show 
the  maximum  response  of  hypothetical  single  degree  of  freedom  resonators  with  different 
natural  frequencies  and  damping  ratios.  In  this  report  response  spectra  were  calculated  for  a 
5%  of  critical  damping  ratio.  The  spectra  are  presented  in  two  plots.  One  plot  compares 
acceleration  response  spectra,  plotted  as  a function  of  period  on  a semi-logarithmic  scale 
which  magnifies  the  high  frequency  (low  period)  range.  A second  plot  compares  relative 
displacements,  pseudo  relative  velocities,  and  pseudo  accelerations  as  a function  of  period  in 
a single  plot  using  four  logarithmic  scales.  This  plot  conveys  more  information  on  the  spectra. 
The  plot  uses  virtual  (pseudo)  values  for  relative  velocities  and  accelerations,  using  the 
following  relationships,  which  are  for  sinusoidal  forcing  functions,  and  produce  quantities 
which,  on  the  average,  tend  to  be  smaller  than  those  actually  calculated  for  the  time  histories: 

PRV=U(  2 a/)  • • • (4) 


PAA  = U(2-w  f)2  • • • (5) 


Where:  PRV  = pseudo  spectral  relative  velocity 

PAA  = pseudo  spectral  absolute  acceleration 
U = calculated  spectral  relative  displacement  amplitude 
f = frequency,  Hz  (reciprocal  of  period  in  plot) 


In  accordance  with  conventions  for  building  design,  periods  rather  than  frequencies  are  shown 
in  the  plots.  In  some  instances  Fourier  spectra  are  shown  in  order  to  provide  more  information 
about  the  frequency  content  of  acceleration  time  histories. 
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4.2 


Discussion  of  Results 


4.2.1  Oakland  Outer  Harbor  Wharf  Motion. 

Some  comparison  of  a predicted  with  an  observed  acceleration  time  history  can  be  obtained 
for  the  Loma  Prieta  strong  motion  record  at  the  Oakland  Outer  Harbor  Wharf.  Calculated 
ground  surface  motions  were  obtained  by  using  the  soil  profile  in  figure  2,  and  the  rock 
outcrop  motion  in  the  East-West  direction  recorded  at  the  Yorba  Buena  strong  motion  station. 

In  this  instance,  a comparison  of  the  recorded  Oakland  Outer  Harbor  Wharf  acceleration  time 
history  with  the  calculated  time  history  of  the  ground  surface  motion  indicated  that  there  was 
a tendency  to  underestimate  the  amplitude  of  the  ground  motion.  Since  some  of  the  variables 
used  to  calculate  the  ground  surface  motion  are  estimated,  a range  of  values  for  these 
estimated  variables  was  used  in  order  to  ascertain  their  effect  on  the  calculated  ground 
surface  motion.  The  estimated  variables  include  the  baserock  stiffness  (for  which  no  data 
were  available),  the  rock  outcrop  motion  used  as  input  for  the  calculations  (the  Yorba  Buena 
and  the  Rincon  Hill  records  were  used),  and  the  effective  strain  used  in  the  equivalent  linear 
model  (a  quantity  which  is  specified). 

Figure  4.1  shows  three  acceleration  time  histories:  the  Yorba  Buena  rock  outcrop  motion,  and 
time  histories  at  the  soil-rock  interface,  calculated  for  base  rock  stiffnesses  associated  with 
a mass  density  of  the  base  rock  of  2.56  Mg/m3  and  shear  wave  propagation  velocities  in  the 
base  rock  of  762  m/s  and  1 ,220  m/s,  respectively.  The  time  histories  at  the  soil/rock  interface 
for  these  two  latter  cases  were  calculated  assuming  that  the  time  history  of  the  incident  wave 
for  these  cases  was  identical  with  the  time  history  of  the  incident  wave  which  produced  the 
Yorba  Buena  rock  outcrop  motion.  It  can  be  seen  from  figure  4.1  that  the  amplitude  of  the 
motion  at  the  soil/rock  interface  increases  with  the  stiffness  of  the  base  rock.  (For  the  case 
of  Vs  (rock)  = 762  m/s  the  maximum  acceleration  at  the  soil/rock  interface  was  56%  of  that 
of  the  outcrop  motion,  and  for  the  case  of  V5(rock)  = 1220  m/s  it  was  65%  of  that  of  the 
outcrop  motion).  An  upper  bound  for  the  motion  at  the  soil/rock  interface,  which  would  occur 
for  an  infinitely  stiff  base  rock,  would  be  a motion  at  the  soil/rock  interface  which  is  identical 
with  the  Yorba  Buena  rock  outcrop  motion.  It  is  assumed  in  this  latter  analysis  that  the 
bedrock  motion  at  Yorba  Buena  was  identical  with  that  at  the  base  of  the  deposit,  except  for 
the  effect  of  confinement.  This  was  not  necessarily  the  case,  and  it  is  not  unreasonable  to 
assume  that  some  of  the  energy  of  the  bedrock  motion  was  lost  as  it  propagated  to  the 
exposed  bedrock  surface,  which  is  located  at  a much  higher  elevation.  Figure  4.2  shows 
three  time  histories:  the  time  history  recorded  at  the  Oakland  Outer  Harbor  Wharf  strong 
motion  station  in  a 305°  direction;  a calculated  ground  motion  using  the  Yorba  Buena  rock 
outcrop  motion  in  an  East-West  direction  at  the  soil/rock  interface  (using  the  unconfined  rock 
motion  at  the  base  of  the  deposit  is  equivalent  to  the  assumption  that  the  rock  is  infinitely 
stiff);  and  a calculated  ground  motion  using  a confined  rock  motion  calculated  from  the  Yorba 
Buena  rock  outcrop  motion  and  assuming  that  the  shear  wave  propagation  velocity  in  the  base 
rock  is  1 ,200  m/s. 

Figures  4.3  and  4.4  compare  response  spectra  calculated  for  5%  damping.  The  calculated 
response  spectra  are  for  the  recorded  ground  motion,  and  for  calculated  ground  motions  using 
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the  Yorba  Buena  rock  outcrop  motion  at  the  soil/rock  interface,  and  confined  baserock 
motions  calculated  for  base  rock  shear  wave  propagation  velocities  of  1200  and  762  m/s. 

Figure  4.3  compares  acceleration  response  spectra,  and  figure  4.4  compares  response  spectra 
for  relative  displacement,  pseudo  relative  velocity,  and  pseudo  acceleration.  It  is  interesting 
to  note,  that  in  terms  of  frequency  response  the  spectra  are  quite  similar,  however  they  differ 
in  amplitude.  It  is  evident  that  the  set  of  assumptions  used  in  the  calculations  tended  to 
underestimate  the  amplitude  of  the  spectral  response  for  periods  less  than  2 s. 

For  the  sake  of  comparison,  calculations  were  also  performed  using  the  Rincon  Hill  outcrop 
motion  record,  which  is  not  as  close  to  the  Oakland  Outer  Harbor  Wharf  site  as  the  Yorba 
Buena  outcrop  motion.  This  is  shown  in  figures  4.5  and  4.6  for  an  assumed  shear  wave 
propagation  velocity  in  the  base  rock  of  762  m/s.  It  can  be  seen  from  the  figure,  that  the 
results  obtained  from  the  Yorba  Buena  and  the  Rincon  Hill  outcrop  motions  are  similar,  and 
that  those  calculated  from  the  Yorba  Buena  motion  more  closely  resemble  the  frequency 
response  obtained  from  the  recorded  ground  motion. 

Another  variable  that  was  investigated  is  the  effective  strain  used  to  calculate  the  ground 
motion,  which  is  specified  as  a percentage  of  the  calculated  maximum  strain  for  each  layer. 
The  effective  strain  is  used  to  select  dynamic  soil  properties  which  in  turn  are  used  to 
calculate  the  ground  motion.  Three  effective  strain  levels  were  checked:  55%  of  maximum, 
which  was  used  in  most  of  the  calculations  in  this  study;  61%  of  maximum  which 
corresponds  to  an  effective  strain  level  of  [(M-l)/10]  ^max,  where  M = Richter  Magnitude  of 
Earthquake  and  7^  is  the  maximum  shear  strain  calculated  for  the  time  history  in  the  center 
of  each  layer,  a strain  level  suggested  in  another  part  of  the  NIST  study  (refer  to  introduction), 
and  35.5%  of  maximum,  which  corresponds  to  the  ratio  of  the  RMS  value  of  earthquake 
accelerations,  calculated  between  the  time  limits  of  0.1  5 and  0.85  in  the  Husid  plot,  to  the 
maximum  acceleration.  [The  Husid  plot  is  a plot  of  equation  6 below  and  is  indicative  of  the 
fraction  of  the  earthquake  energy  released  as  a function  of  time  (Idriss,  1 978,  Husid,  1 968)]. 


J A2(t)dt 

0 

f A2(t)dt 

0 


(6) 


where:  A = Acceleration  at  time  t 

tf  = time  at  end  of  record 

The  effect  of  the  level  of  the  assumed  effective  strain  is  shown  in  figures  4.7  and  4.8.  It  is 
evident  that,  within  the  range  examined,  the  level  of  effective  strain  did  not  significantly  affect 
the  results  of  the  calculations.  For  assumed  effective  strain  levels  of  35.5%,  55%,  and  61  %, 
maximum  calculated  strains  were  0.17%,  0.23%,  and  0.25%,  respectively,  in  the  center  of 
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the  most  strained  layer,  calculated  dominant  periods  were  1.69s,  1.75s,  and  1.77s, 
respectively,  and  maximum  accelerations  at  the  surface  of  the  deposit  were  0.21 8 G,  0.204G, 
and  0.1 98G,  respectively.  It  can  be  seen  from  figure  1 , that  for  the  strain  levels  in  question 
a change  in  strain  of  0.08  percent  will  not  substantially  affect  the  dynamic  soil  properties. 

Similarly,  use  of  dynamic  soil  properties  recommended  by  Sun,  1 988,  and  Vucetic  and  Dobry, 
1991,  had  little  effect  on  the  amplitude  of  the  calculated  ground  motion. 

In  summary,  calculations  using  equivalent  linear  response  analysis  and  incorporating  the 
dynamic  soil  properties  shown  in  figure  2.1 , preliminary  soil  deposit  data  from  the  site  of  the 
Oakland  Outer  Harbor  Wharf,  and  base  rock  motions  equivalent  to  the  Yorba  Buena  and 
Rincon  Hill  rock  outcrop  motions  recorded  in  the  Loma  Prieta  Earthquake,  produced  ground 
motions  at  the  surface  of  the  deposit  which  tended  to  be  smaller  in  amplitude  but  reasonably 
similar  in  frequency  content  when  compared  with  the  strong  motion  time  history  recorded 
during  the  Loma  Prieta  earthquake.  An  increase  in  the  assumed  base  rock  stiffness,  for  which 
no  data  were  available,  tended  to  reduce,  but  not  eliminate,  the  difference.  However,  the 
difference  is  eliminated  when  the  base  rock  is  assumed  to  be  rigid.  This  latter  assumption 
cannot  be  justified. 

It  is  suggested  that  the  amplitude  of  the  incident  shear  waves  at  the  base  of  the  deposit  may 
have  been  greater  than  that  recorded  at  the  Yorba  Buena  and  Rincon  Hill  strong  motion 
stations.  Another  possible  reason  for  the  difference  between  calculated  and  observed  ground 
motion  amplitudes  could  be  that  the  damping  ratios  used  in  the  analysis,  which  are  based  on 
laboratory  test  data,  rather  than  in  situ  measurements,  significantly  overestimate  the  actual 
damping  by  the  deeper  soil  deposits.  Oversimplifications,  inherent  in  the  SHAKE  model,  may 
also  have  contributed  to  the  difference. 
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Acceleration,  G's  Acceleration,  G's  Acceleration,  G's 


Yerba  Buena  Rock  Outcrop  Motion 


T~p~r~i  i ttt  i v|  i i i i i i i i i | i i i n~m  ;-j 

10  20  30  40 


Time,  8 


Calculated  Motion  at  Soil/Rock  Interface  Assuming 
a Shearwave  Velocity  of  1220  m/s  in  Underlying  Rock 


Calculated  Motion  at  Soil/Rock  Interface  Assuming 
a Shearwave  Velocity  of  762  m/s  in  Underlying  Rock 


Time,  s 


Figure  4. 1 Effect  of  Base  Rock  Stiffness  on  the  Base  Rock  Motion  amplitude  at  the  Oakland 
Outer  Harbor  Wharf  Strong  Motion  Site,  Calculated  for  the  Yorba  Buena  Rock 
Outcrop  motion. 
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Figure  4.2: 


Calculated  Ground  Motion  Using  the  Yerba 
Buena  Record  as  Soil/Rock  Interface  Motion 


Calculated  Ground  Surface  Motion  Using  the  Yerba 
Buena  Record  as  an  Outcrop  Motion  and  Assuming  a 
Bedrock  Shearwave  Propagation  Velocity  of  1220  m/s 


Recorded  and  Calculated  Acceleration  Time  History  at  the  Oakland 
Outer  Harbor  Wharf  Strong  Motion  Station. 
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Acceleration,  G's 


Figure  4.3:  Comparison  of  Acceleration  Response  Spectra  for  5%  Damping  for  the 
Recorded  Ground  Motion  and  for  Calculated  Ground  Motions  at  the  Oakland 
Outer  Harbor  Wharf  Strong  Motion  Station. 
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Pseudo  Relative  Velocity,  mm/s 


Recorded  Oakland  Wharf  Ground  Motion 
Yorba  Buena  Rock  Outcrop  Motion 


Calculated  Ground  Surface  Motions 


— — — Yorba  Buena  Outcrop  Motion  at  Soil/rock  Interface 
— — Confined  Rock  Motion,  Vs  frock}  = 1,200  m/s 
.... — ..  Confined  Rock  Motion,  Vs  (rock)  = 760  m/s 
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Figure  4.4:  Comparison  of  Displacement,  Pseudo  Relative  Velocity,  and  Pseudo 
Acceleration  Response  Spectra  for  5%  Damping  for  the  Recorded 
Ground  Motion  and  for  Calculated  Ground  Motions  at  the  Oakland 
Outer  Harbor  Wharf  Strong  Motion  Station. 
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Acceleration,  G's 


1 1 Oakland  Wharf  Ground  Motion 

—  Rincon  Hill  Bedrock  Motion 

—  Yorba  Buena  Bedrock  Motion 

Calculated  Ground  Surface  Motions 

— — «•  Case  14,  Yorba  Buena  Outcrop  Motion 
Case  14,  Rincon  Hill  Outcrop  Motion 


Figure  4.5:  Comparison  of  Acceleration  Response  Spectra  for  5%  Damping  for  the  Recorded 
Oakland  Outer  Harbor  Wharf  Ground  Motion  and  for  Ground  Motions  Calculated 
Using  the  Yorba  Buena  and  the  Rincon  Hill  Rock  Outcrop  Motions. 
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Pseudo  Relative  Velocity,  mm/s 


■ Recorded  Oakland  Wharf  Ground  Motion 

Rincon  Hill  Outcrop  Motion 

— — - Yorba  Buena  Outcrop  Motion 

Calculated  Ground  Surface  Motions 

" - - Case  14,  Yorbo  Buena  Outcrop  Motion 
— — — Case  14,  Rincon  Hill  Outcrop  Motion 
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Figure  4.6:  Comparison  of  Response  Spectra  for  5%  Damping  for  the  Recorded 
Oakland  Outer  Harbor  Wharf  Ground  Motion  and  for  Ground  Motions 
Calculated  Using  the  Yorba  Buena  and  the  Rincon  Hill  Rock  Outcrop  Motions. 
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Acceleration,  G*s 


Recorded  Oakland  Wharf  Ground  Motion 
Colculoted  Ground  Surface  Motions 

----  Calculated  Motion,  Effective  Strain  35.5*  of  Maximum 

Calculated  Motion,  Effective  Strain  55*  of  Maximum 

Calculated  Motion,  Effective  Strain  61*  of  Maximum 


Figure  4.7:  Effect  of  Assumed  Effective  Strain  Level  on  the  Calculated  Ground  Motion 
Acceleration  Response  Spectra. 
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Figure  4.8:  Effect  of  Assumed  Effective  Strain  Level  on  Response  Spectra  for  the  Calculated 
Ground  Motion. 
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4.2.2  Profiles  with  New  Bay  Mud  (Cases  1 through  4) 

4.2.2. 1 Profiles  Studied 

The  profiles  are  shown  in  Table  1.  Profile  1 has  3m  (10  ft)  of  fill,  3 m of  New  Bay  Mud,  3 m 
of  alluvium,  and  3 m of  Old  Bay  Mud  over  a halfspace  of  bedrock.  The  sequence  of  these 
layers  is  in  the  same  order  as  it  is  encountered  in  many  locations  in  the  San  Francisco  Bay 
area.  In  profiles  2,  3,  and  4,  the  depth  of  the  New  Bay  mud  is  increased  to  6,  1 5,  and  30  m, 
respectively.  The  intent  of  these  hypothetical  profiles  is  to  study  the  effect  of  increasing  the 
depth  of  the  New  Bay  Mud.  In  all  the  cases  studied  the  groundwater  level  is  at  the  bottom 
of  the  fill,  and  the  dynamic  shear  modulus  of  the  bedrock  at  small  strains  is  assumed  to  be 
2,320  MPa. 

4. 2. 2. 2 Effect  of  New  Bay  Mud  depth 

Response  Spectra  for  a 5%  critical  damping  ratio  for  ground  motions,  calculated  for  cases  1 , 
3,  and  4 when  subjected  to  a base  rock  excitation  represented  by  the  Rincon  Hill  outcrop 
motion  in  the  East-West  direction,  are  shown  in  figures  4.9  and  4.10.  The  amplification  at 
very  low  periods  shown  in  figure  4.9  is  a function  of  the  maximum  acceleration  in  the  time 
history  of  the  calculated  ground  motions  (at  very  low  periods  the  spectral  acceleration 
response  equals  the  maximum  acceleration).  It  can  be  seen  from  figure  4.9,  that  the  maximum 
acceleration  of  the  time  history  was  amplified  by  all  the  soil  profiles.  The  maximum 
acceleration  of  the  Rincon  Hill  record  is  0.09G.  This  compares  with  A ^ values  of  0.1 7G, 
0.1 9G,  and  0.21  G for  cases  1, 3,  and  4,  respectively,  which  represent  amplifications  with 
respect  to  the  rock  outcrop  motions  of  1.9,  2.1,  and  2.3,  respectively.  This  compares  with 
average  amplification  ratios  of  2.6  observed  for  Bay  Mud  sites  in  the  Loma  Prieta  Earthquake 
(EERI,  1989).  The  ratios  between  the  maximum  acceleration  of  the  base  rock  motions  and 
that  of  the  rock  outcrop  motion  is  0.89,  0.86,  and  0.83,  respectively,  for  cases  1 , 3,  and  4. 

As  can  be  seen  from  figures  4.9  and  4.10,  the  trend  of  the  effect  of  the  depth  of  the  New 
Bay  Mud  is  to  amplify  the  response  in  the  vicinity  of  the  characteristic  period  of  the  deposit. 
Thus  profiles  3 and  4 result  in  a high  response  for  periods  ranging  from  0.5  to  2 seconds. 
This  is  the  range  of  periods  which  will  affect  most  low  to  medium  rise  buildings.  Case  2 is  an 
interesting  exception  to  the  otherwise  consistent  trend  of  the  calculated  results.  Figure  4.1 1 
shows  the  data  from  figure  4.9,  with  the  acceleration  response  spectrum  for  Case  2 
superimposed  and  highlighted.  The  maximum  calculated  ground  acceleration  for  Case  2 is 
only  0.1 3G,  which  represents  an  amplification  of  A ^ of  1.4  with  respect  to  the  Rincon  Hill 
outcrop  motion.  Note  that  the  characteristic  period  for  Profile  2 is  0.39s.  This  coincides  with 
a low  spectral  response  for  the  Rincon  Hill  motion. 

The  cyclic  strains  induced  in  the  various  soil  layers  by  the  Rincon  Hill  rock  outcrop  motion  are 
relatively  small.  Maximum  shear  strains  range  from  0.07%  in  Profile  2 to  0.14%  in  Profile  4. 
It  is  suggested  that  the  SHAKE  program  could  be  used  up  to  maximum  cyclic  strains  of  the 
order  of  1 %,  provided  the  soils  are  not  susceptible  to  liquefaction.  Thus  another  strong  motion 
record  was  selected  to  investigate  the  response  in  the  vicinity  of  the  suggested  upper  limit 
of  cyclic  strain.  The  earthquake  motion  selected  was  the  rock  outcrop  motion  in  the  North- 
South  direction  recorded  for  the  Loma  Prieta  earthquake  by  the  strong  motion  station  in  Santa 
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Cruz.  This  latter  strong  motion  record,  which  is  for  an  outcropping  bedrock  site,  has  a 
maximum  acceleration  of  0.441 3G  and  a frequency  content  which  differs  from  that  of  the 
Rincon  Hill  earthquake.  The  difference  in  frequency  content  is  illustrated  in  figures  4.12  and 
4.13.  Figure  4.12  compares  Fourier  spectra,  and  figure  4.13  normalized  response  spectra.  It 
is  thought  that,  at  least  in  part,  the  difference  in  the  frequency  content  of  these  ground 
motions  is  attributable  to  their  distance  from  the  source  of  the  ground  motion.  Figure  4.13 
distorts  the  relative  scale  of  the  response  spectra.  In  figure  4.14  the  same  spectra  are  plotted 
to  scale. 

Figures  4.15  and  4.16  show  a comparison  between  acceleration  response  spectra  for  the 
recorded  Santa  Cruz  motion  and  those  calculated  at  the  ground  surface  for  profiles  1 , 3,  and 
4.  It  can  be  seen  that  while  the  response  in  the  low  period  (high  frequency)  range  is 
substantially  amplified  by  profile  1 , there  is  no  amplification  for  periods  smaller  than  0.2s  by 
profiles  3 and  4.  The  difference  between  the  responses  to  the  low  amplitude  Rincon  Hill 
motion  and  the  high  amplitude  Santa  Cruz  motion  is  best  illustrated  by  the  calculated 
amplification  spectra  for  these  motions.  This  is  shown  in  figure  4.1 7.  The  amplification  factor 
plotted  in  the  figure  is  the  ratio  between  the  amplitudes  of  the  ground  surface  motion  and  that 
of  the  base  rock  motion,  which  is  a function  of  frequency.  In  this  case  the  maximum 
acceleration  of  the  base  rock  motion  was  approximately  85%  of  that  of  the  rock  outcrop 
motion.  Note  that  for  the  Santa  Cruz  earthquake  the  amplification  factor  dropped  below  unity 
for  frequencies  in  excess  of  5 Hz,  while  for  the  Rincon  Hill  it  dropped  below  unity  for 
frequencies  in  excess  of  10  Hz.  In  this  instance  Profile  2 followed  the  trend  of  the  other 
profiles.  This  is  illustrated  in  figure  4.1 8,  where  the  acceleration  response  spectrum  for  Case 
2 is  superimposed  on  figure  4.15  and  highlighted. 

It  is  also  of  interest  to  compare  the  response  spectra  for  case  4,  calculated  by  the  Rincon  Hill 
and  Santa  Cruz  motions.  This  is  done  in  figures  4. 1 9 and  4.20.  Note  that  for  periods  in  excess 
of  1 .3s  the  amplitudes  of  the  response  spectra  for  these  two  bedrock  motions  are  about 
equal.  The  implication  of  this  observation  would  be  that  in  deep  soft  soils  the  response  of 
taller  buildings  would  not  be  very  sensitive  to  the  distance  of  the  source  of  the  ground  motion. 
However  this  phenomenon  is  at  least  in  part  attributable  to  the  difference  in  the  frequency 
contents  between  the  two  bedrock  motions.  While  this  difference  appears  to  be  logical,  it  has 
not  been  demonstrated  at  this  time  that  it  represents  a trend  which  is  typical  for  most 
earthquake  data. 

4. 2. 2. 3 Effect  of  Layer  Thickness  used  in  Analysis 

As  previously  noted,  the  revised  SHAKE  program  developed  in  this  project  enables  us  to  use 
50  soil  layers  (including  the  halfspace).  This  enables  us  to  accommodate  deep  and  complex 
soil  profiles,  such  as  that  associated  with  the  Oakland  Outer  Harbor  Wharf,  which  was 
previously  discussed.  It  was  thought  that  a larger  number  of  layers  would  also  be  useful  in 
studying  the  high  frequency  range  of  the  spectrum.  It  was  reasoned  that  layer  thickness 
should  not  exceed  1/4  of  the  wavelength  associated  with  the  highest  frequency  considered 
in  the  analysis.  In  this  study,  all  frequencies  above  25  Hz  were  filtered  out.  Thus, 
theoretically,  layer  depths  should  correspond  to  the  wavelengths  associated  with  this 
frequency.  In  practice,  however,  for  the  time  histories  used  in  this  analysis,  very  small 
acceleration  amplitudes  were  associated  with  frequencies  above  15  Hz.  Thus  it  would  be 
difficult  to  determine  the  effect  of  layer  thickness  on  the  basis  of  calculated  time  histories  or 
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response  spectra  (on  the  high  frequency  end  of  the  response  spectra  the  acceleration 
response  amplitude  simply  equals  the  maximum  acceleration  in  the  record  and  is  not  sensitive 
to  the  high  frequency  components  of  the  input  motion).  A more  direct  way  to  determine  the 
effect  of  layer  thickness  is  therefore  via  the  amplification  spectra. 

Figures  4.21  through  4.23  illustrate  the  effect  of  varying  the  layer  thickness.  Figure  4.21 
represents  a shallow  soil  profile  with  relatively  few  layers.  In  this  instance  the  effect  of 
increasing  the  number  of  layers  was  substantial.  It  should  be  noted  however,  that  the  effect 
on  the  amplification  spectrum,  which  is  confined  to  frequencies  over  7Hz,  did  not  show  up 
in  the  calculated  response  spectra.  In  figures  4.22  and  4.23,  which  represent  deep  soil 
profiles,  the  effect  of  changing  the  number  of  layers  is  insignificant.  It  was  therefore 
concluded  that  for  deeper  soil  profiles  it  was  not  necessary  to  keep  the  layer  thickness  below 
1/4  wave  length. 

4.2.3  Effect  of  Alluvium  (Cases  5 to  7) 

Alluvial  profiles  are  explored  in  cases  5 to  7.  Response  spectra  for  cases  5 to  7 are  shown  in 
figures  4.24  and  4.25  for  the  Rincon  Hill  rock  outcrop  motion  in  the  East  West  direction,  and 
in  figures  4.26  and  4.27  for  the  Santa  Cruz  rock  outcrop  motion  in  the  North-South  direction. 
The  spike  in  the  response  spectrum  for  profile  6 in  figure  4.24  is  attributable  to  the 
characteristic  period  of  profile  6 of  0.27s  which  coincides  with  a peak  in  the  response 
spectrum  for  the  Rincon  Hill  outcrop  motion.  Amplification  of  the  calculated  ground  surface 
motions  for  Profiles  5 to  7 with  respect  to  the  maximum  acceleration  of  the  Rincon  Hill 
outcrop  motion  range  from  1 .45  to  1 .8.  This  compares  with  an  average  amplification  factor 
of  1 .8  for  alluvial  sites  in  the  San  Francisco  Bay  area  observed  in  the  Loma  Prieta  earthquake 
(EERI,  1989). 

Note  that  for  the  Santa  Cruz  outcrop  motion  the  response  spectra  for  profiles  6 and  7 are 
virtually  identical  between  the  periods  of  0.5  and  1 .2s,  which  are  of  interest  with  respect  to 
building  response.  There  was  substantial  amplification  of  the  Santa  Cruz  outcrop  motion  for 
periods  below  1 .2s.  For  both,  the  Santa  Cruz  and  the  Rincon  Hill  motion  there  was  no 
amplification  for  periods  exceeding  1 .2s.  Amplification  spectra  for  case  7 for  the  Rincon  Hill 
and  Santa  Cruz  motions  are  compared  in  figure  4.28. 

Response  spectra  for  cases  4 and  7 are  compared  in  figures  4.29  through  4.32.  Note  that  for 
the  Rincon  Hill  outcrop  motion  the  low  period  range  of  the  spectrum  is  about  equally  amplified 
by  profiles  4 and  7,  and  that  the  Santa  Cruz  motion  is  substantially  amplified  in  the  low  period 
range  by  profile  7. 

4.2.4  Composite  Profiles 

4.2.4. 1 New  Bay  Mud  Over  Alluvium 

The  effect  of  alluvial  deposits  underlying  New  Bay  mud  was  explored  by  profiles  8 and  9.  The 
results  are  illustrated  in  figures  4.33  through  4.38.  In  figures  4.33  and  4.34  cases  7,  9,  and 
4 are  compared  for  the  Rincon  Hill  outcrop  motion.  Note  that  there  is  no  significant  difference 
between  cases  9 and  4.  Thus  the  addition  of  a 24  m thick  layer  of  alluvium  at  the  base  of 
the  30  m layer  of  New  Bay  mud  did  not  significantly  alter  the  calculated  effect  of  the  ground 
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surface  motion  on  structures.  A similar  conclusion  can  be  drawn  from  figures  4.34  and  4.35, 
which  compare  the  response  of  profiles  7,  9,  and  4 for  the  Santa  Cruz  outcrop  motion. 
Amplification  spectra  for  cases  9 and  4 are  compared  in  figures  4.37  and  4.38  for  the  Rincon 
Hill  and  Santa  Cruz  outcrop  motions,  respectively. 

4. 2. 4. 2 Old  Bay  Mud,  Alluvium  and  New  Bay  Mud 

The  combination  of  Old  Bay  mud  (or  comparably  stiffer  soils),  alluvium,  and  New  Bay  mud, 
which  is  common  in  the  Bay  area,  is  explored  in  profile  1 3,  which  has  1 50  m of  Old  Bay  mud, 
6 m of  alluvium,  6 m of  New  Bay  mud,  and  3 m of  fill.  This  profile  is  thought  to  be  somewhat 
similar  to  that  of  the  site  of  the  1-880  collapse  in  Oakland.  In  figures  4.39  and  4.40  response 
spectra  for  calculated  ground  motions  for  Case  1 3,  using  a maximum  base  rock  shear  wave 
velocity  of  1 ,220  m/s,  are  compared  with  response  spectra  for  the  recorded  ground  motions 
in  the  East  West  direction  at  Emeryville  and  Oakland  Outer  Harbor  Wharf.  The  Emeryville 
ground  motion,  which  has  a response  peak  at  the  1.0-s  period,  could  have  been  influenced 
by  an  adjacent  building.  If  this  was  the  case,  it  would  not  be  a true  free  field  motion.  This  is 
illustrated  by  the  regularity  of  the  calculated  velocity  and  displacement  time  histories  which 
are  shown  in  figure  4.41 . The  correlation  between  these  calculated  ground  motions  and  the 
recorded  ground  motion  at  Oakland  Outer  Harbor  Wharf  is  similar  to  that  with  the  ground 
motion  which  was  calculated  from  actual  subsurface  data.  In  figure  4.42,  the  calculated 
ground  motion  at  Oakland  Outer  Harbor  Wharf  is  compared  with  the  ground  motions 
calculated  for  case  13.  Note  that  for  periods  greater  than  0.6s  there  is  no  significant 
difference  between  the  three  calculated  ground  motions.  Thus  the  hypothetical  soil  profiles 
used  in  this  study  seem  to  produce  results  which  are  similar  to  those  obtained  from  actual  soil 
data.  However,  in  both  instances  the  amplitude  of  the  ground  motion  was  underestimated. 

4.2.5  Comparison  of  New  Bay  Mud,  Alluvium,  and  Old  Bay  Mud 

In  profiles  4,  11,  and  12,  the  thickness  of  New  Bay  mud,  alluvium,  and  Old  Bay  mud  is 
increased  to  30  m,  while  the  thickness  of  the  other  layers  in  profile  1 is  kept  at  3 m.  Cases 
4,  1 1,  and  12,  calculated  for  the  Rincon  Hill  and  Santa  Cruz  outcrop  motions  are  compared 
in  figures  4.43  to  4.46.  Predictably,  the  results  for  the  alluvium  and  Old  Bay  mud  are  quite 
similar. 
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Figure  4.9:  Comparison  of  Acceleration  Response  Spectra  for  Cases  1, 3,  and  4 - Rincon 
Hill  Outcrop  Motion. 
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Case  4,  30  m New  Bay  Mud 
Case  3,  15  m New  Bay  Mud 
Case  1,  3 m New  Bay  Mud 
Rincon  Hill  Outcrop  Motion 
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Figure  4.10:  Comparison  of  Response  Spectra  for  Cases  1,  3,  and  4 - Rincon  Hill 
Outcrop  Motion. 
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Acceleration,  G's 


Figure  4.11:  Comparison  of  Acceleration  Response  Spectrum  for  Case  2 
those  for  Cases  1, 3,  and  4 - Rincon  Hill  Outcrop  Motion. 
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Figure  4.12:  Fourier  Spectra  for  the  Rincon  Hill  and  Santa  Cruz  Rock  Outcrop  Motions. 
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Figure  4.13: 


Comparison  of  Normalized  Acceleration  Response  Spectra  for  the 
Rincon  Hill  and  Santa  Cruz  Rock  Outcrop  Motions. 
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Figure  4.14: 


Comparison  of  Acceleration  Response  Spectra  for  the  Rincon  Hill  and  Santa 
Cruz  Rock  Outcrop  Motions. 
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Acceleration,  G's 


Figure  4.15:  Comparison  of  Acceleration  Response  Spectra  for  Cases  1,  3, 

and  4 - Santa  Cruz  Rock  Outcrop  Motion. 
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Figure  4.16:  Comparison  of  Response  Spectra  for  Cases  1 and  4 - Santa  Cruz 
Rock  Outcrop  Motion. 
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Figure  4.17:  Amplification  Spectra  for  Case  4. 
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Figure  4.1 8:  Comparison  of  Acceleration  Response  Spectra  for  Cases  1 , 2,  3,  and  4 - Santa 
Cruz  Rock  Outcrop  Motion. 
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Cas«  4.  Sonto  Cruz  Outcrop  Motion 
— — — Cos*  4,  Rincon  Hill  Outcrop  Motion 


Figure  4.19:  Comparison  of  Acceleration  Response  Spectra  for  Profile  4,  Rincon  Hill  and 
Santa  Cruz  Outcrop  Motions. 
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Figure  4.20:  Comparison  of  Response  Spectra  for  Profile  4,  Rincon  Hill  and 
Santa  Cruz  Outcrop  Motions. 
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Figure  4.21:  Amplification  Spectra  for  Case  1,  Calculated  for  9 and  1 1 Layers. 
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Figure  4.22:  Amplification  Spectra  for  Case  3,  Calculated  for  10  and  27  Layers. 


42 


Frequency,  Hz 


Figure  4.23:  Amplification  Spectra  for  Case  4,  Calculated  for  25  and  45  Layers. 
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Figure  4.24:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  5 to  7,  Calculated  for 
the  Rincon  Hill  Outcrop  Motion. 
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Figure  4.25:  Comparison  of  Response  Spectra  for  Profiles  5 to  7,  Calculated  for  the  Rincon 
Hill  Outcrop  Motion. 
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Figure  4.26:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  5 to  7,  Calculated  for 
the  Santa  Cruz  Outcrop  Motion. 
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Figure  4.27:  Comparison  of  Response  Spectra  for  Profiles  5 to  7,  Calculated  for  the  Santa 
Cruz  Outcrop  Motion. 
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Figure  4.28:  Comparison  of  Amplification  Spectra  for  Case  7,  Calculated  for  the  Rincon  Hill 
and  Santa  Cruz  outcrop  motions. 
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Figure  4.29:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  4 and  7,  Calculated 
for  the  Rincon  Hill  Outcrop  Motion. 
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Figure  4.30:  Comparison  of  Response  Spectra  for  Profiles  4 and  7,  Calculated  for  the  Rincon 
Hill  Outcrop  Motion. 
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Figure  4.31 : Comparison  of  Acceleration  Response  Spectra  for  Profiles  4 and  7,  Calculated 
for  the  Santa  Cruz  Outcrop  Motion. 
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Figure  4.32:  Comparison  of  Response  Spectra  for  Profiles  4 and  7,  Calculated  for  the  Santa 
Cruz  Outcrop  Motion. 
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Figure  4.33:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  4,  7, 

and  9,  Calculated  for  the  Rincon  Hill  Outcrop  Motion. 
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Figure  4.34:  Comparison  of  Response  Spectra  for  Profiles  4,  7,  and  9,  Calculated  for  the 
Rincon  Hill  Outcrop  Motion. 
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Figure  4.35:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  4,  7, 

and  9,  Calculated  for  the  Santa  Cruz  Outcrop  Motion. 
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Figure  4.36:  Comparison  of  Response  Spectra  for  Profiles  4,  9,  and  7,  Calculated  for  the 
Santa  Cruz  Outcrop  Motion. 
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Figure  4.37:  Comparison  of  Amplification  Spectra  for  Cases  4 and  9,  Rincon  Hill  Outcrop 
Motion. 
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Figure  4.38:  Comparison  of  Amplification  Spectra  for  Cases  4 and  9,  Santa  Cruz  outcrop 
motion. 
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Figure  4.39:  Comparison  of  Acceleration  Response  Spectra  for  Case  13  with  the  Recorded 
Ground  Motions  at  Oakland  Outer  Harbor  Wharf  and  Emeryville. 
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Figure  4.40:  Comparison  of  Response  Spectra  for  Case  13  with  the  Recorded  Ground 
Motions  at  Oakland  Outer  Harbor  Wharf  and  Emeryville. 
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Figure  4.41:  Acceleration,  Velocity,  and  Displacement  Time  Histories  for  the  Recorded 
Emeryville  Ground  Motion. 
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Figure  4.42:  Comparison  of  Calculated  Ground  Motions  for  Case  13  and  Oakland  Outer 
Harbor  Wharf. 
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Figure  4.43:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  4,  1 1 
Calculated  for  the  Rincon  Hill  Outcrop  Motion. 
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Figure  4.44:  Comparison  of  Response  Spectra  for  Profiles  4,  1 1 , and  1 2,  Calculated  for  the 
Rincon  Hill  Outcrop  Motion. 
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Figure  4.45:  Comparison  of  Acceleration  Response  Spectra  for  Profiles  4,  11,  and  12 
Calculated  for  the  Santa  Cruz  Outcrop  Motion. 
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Figure  4.46:  Comparison  of  Response  Spectra  for  Profiles  4,  1 1 , and  1 2,  Calculated  for  the 
Santa  Cruz  Outcrop  Motion. 
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5. 


COMPARISON  OF  CALCULATED  RESPONSE  SPECTRA  WITH  NEHRP  DESIGN 
SPECTRA 


5.1  Design  Spectra  Considered 

Elastic  design  spectra  are  used  to  specify  the  base  shear  of  buildings  as  a function  of  their 
fundamental  period.  Several  different  design  spectra  are  presently  used.  These  include  among 
others  the  SEAOC  spectra  (SEAOC,  1990),  the  UBC  spectra  (ICBO,  1988)  and  the  NEHRP 
spectra  (NEHRP,  1988).  These  spectra  are  not  identical.  However,  for  the  sake  of  simplicity, 
this  discussion  is  confined  to  the  NEHRP  recommended  spectra. 

Figures  5.1  and  5.2  show  the  elastic  design  spectra  recommended  in  the  1 988  edition  of  the 
NEHRP  recommended  provisions.  Two  mapped  parameters  are  used  to  specify  these  spectra 
for  various  regions  in  the  U.S.:  Aa,  the  seismic  coefficient  representing  the  effective  peak 
acceleration;  and  Av,  the  seismic  coefficient  representing  the  effective  peak  velocity-related 
acceleration.  Aa  is  used  to  define  the  upper,  cutoff  value  of  the  spectra  (the  horizontal  portion 
of  the  spectra  in  figure  5.1).  Av  is  used  to  define  the  curved  portion  of  the  spectra  in  figure 
5.1.  The  rationale  for  specifying  Aa  and  Av  independently  is  that  the  longer-period  components 
of  the  earthquake  motion,  which  determine  the  part  of  the  design  spectrum  defined  by  Av,  will 
decay  more  slowly  as  the  earthquake  motion  propagates  from  the  source.  Thus  the  contours 
for  Av  differ  from  those  for  Aa.  The  spectra  are  defined  so  that  the  peak  spectral  acceleration 
for  rock  outcrop  sites  corresponds  to  2.5v4a,  and  Aa  represents  the  expected  "effective"  peak 
acceleration  (which  corresponds  to,  but  does  not  necessary  equal,  the  expected  peak 
acceleration  of  the  ground  motion).  According  to  the  commentary  (Vol.  2)  of  NEHRP,  1988, 
spectral  accelerations  for  velocity  were  chosen  so  that  V = 0.762-Av,  where  V = effective 
peak  velocity  in  m/s  and  Av  is  in  G's.  V corresponds  to,  but  does  not  necessarily  equal,  the 
expected  peak  velocity  of  the  ground  motion.  It  is  important  to  note  that,  since  the  publication 
of  NEHRP  1988,  it  has  been  decided  by  ballot  of  the  appropriate  NEHRP  committee  that  the 
maximum  acceleration  for  soil  profile  types  S3  and  S4  should  be  raised  from  2.0-Aa  to  2.5-Aa. 

Spectra  are  defined  for  four  soil  profile  types:  Profile  Type  SI  is  defined  as  rock  (Vs  > 762 
m/s)  and  "stiff"  soil  conditions  (stable  deposits  of  sand,  gravel,  or  stiff  clay  less  than  61  m 
deep);  Profile  Type  S2  is  defined  as  stable  deposits  of  sand,  gravel,  or  stiff  clays  deeper  than 
61  m;  Profile  Type  S3  is  defined  as  9 m or  more  of  soft  to  medium  stiff  clays  with  or  without 
intervening  layers  of  cohesionless  soils;  and  Profile  Type  S4  is  defined  as  more  than  21  m of 
soft  clays  or  silts  (V5  < 122m/s).  There  are  gaps  between  the  defined  soil  profiles.  As  a 
consequence,  many  potential  soil  profiles  do  not  fit  any  of  the  four  profile  types  and  their 
interpretation  is  up  to  the  designer. 

For  the  San  Francisco  Region,  Aa  = Av  = 0.4  G,  and  the  equations  defining  the  elastic  design 
spectrum  for  profile  types  1 and  2 become  (NEHRP,  1988): 


0.485 

0.67 
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Where: 


Cs=  spectral  acceleration 
S = soil  profile  coefficient 
T = period 


S = 1 , 1 .2,  1 .5,  and  2 for  profile  type  SI , S2,  S3,  and  S4,  respectively.  For  profile  types  3 
and  4,  Cs  < 0.8,  rather  than  1 (in  the  published  version  of  NEHRP  '88  - a provision  which  has 
now  been  dropped). 

5.2  Comparison  of  the  NEHRP  design  Spectra  with  Response  Spectra  for  Calculated  and 
Observed  Ground  Motions 

Response  spectra  represent  the  peak  response  of  resonators  with  different  fundamental 
frequencies  and  a given  damping  ratio  to  an  earthquake  ground  motion,  while  design  spectra 
are  used  to  specify  the  design  base  shear  of  structures  as  a function  of  their  dynamic 
characteristics.  Thus  there  is  not  necessarily  a 1 to  1 correspondence  between  calculated 
response  spectra  and  specified  design  spectra.  Nevertheless,  one  should  expect  some 
similarity  between  the  specified  design  spectra  and  the  response  spectra  for  observed  or 
anticipated  ground  motions. 

For  the  sake  of  comparison  of  design  spectra  with  response  spectra  for  recorded  or  calculated 
ground  motions  in  the  San  Francisco  Bay  area,  response  spectra  were  calculated  for  profiles 
S31,  S32,  S41,  S42,  S43,  and  S44  (refer  to  Table  1,  p.  9),  using  the  Rincon  Hill  and  Santa 
Cruz  rock  outcrop  motions  as  input.  These  latter  profiles  are  thought  to  represent  a range  of 
conditions  representative  of  the  S3  and  S4  profiles.  The  Oakland  profile  shown  in  figure  2 in 
this  study  does  not  quite  qualify  for  profile  type  S3  conditions,  but  there  are  27  ft  (8.2  m)  of 
soft  deposits,  and  it  therefore  comes  close  to  the  definition  of  a type  S3  profile.  Another  far 
source  time  history  which  is  of  interest  was  recorded  in  a USGS  array  in  the  vicinity  of 
Redwood  City,  CA,  named  Apeel  2,  approximately  47  km  from  the  source  of  the  Loma  Prieta 
earthquake.  Available  preliminary  information  on  the  subsurface  condition  in  Apeel  2 to  a 
depth  of  30  m (USGS,  1991)  is  shown  in  figure  3.2.  On  the  basis  of  adjacent  information 
from  Foster  City,  CA,  it  is  deduced  that  the  profile  shown  in  figure  3.2  is  probably  on  top  of 
a deep  deposit  of  Old  Bay  mud.  The  Apeel  2 profile  has  a 9 m layer  of  very  soft  New  Bay  mud 
and  thus  meets  the  definition  of  a type  S3  soil  profile.  The  response  spectra  for  the  Apeel  2 
time  history  in  both,  the  North  South  and  the  East  West  direction  exhibit  a very  high  peak 
between  the  1 .Os  and  1.1s  periods.  Profiles  4 and  9 are  considered  type  S4,  and  profiles  7, 
1 1 and  1 2 type  S2.  The  outcrop  motions  in  Santa  Cruz,  Rincon  Hill,  and  Yorba  Buena  are  for 
Type  SI  profiles. 

In  figure  5.3  the  design  spectra  are  compared  with  the  response  spectra  (which  should  be 
roughly  equivalent  to  the  design  spectra)  for  recorded  and  calculated  far  source  ground 
motions  which  form  a critical  envelope.  These  include  the  observed  Oakland  Outer  Harbor 
Wharf  and  Apeel  2 motions,  and  the  calculated  ground  motion  for  profile  4.  The  spectrum  for 
profile  7 is  also  shown  for  comparison.  It  can  be  seen  from  figure  5.3  that  the  response 
spectra  for  the  observed  and  calculated  far-source  ground  motions  in  the  Loma  Prieta 
earthquake  do  not  all  fall  within  the  design  spectra  for  the  San  Francisco  region,  even  if  the 
design  spectra  for  all  the  soil  profile  types  are  extended  to  the  0.8  G maximum  acceleration. 
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Of  particular  concern  is  the  Apeel  2 ground  motion,  which  raises  serious  questions  about  the 
adequacy  of  the  spectrum  for  profile  type  S3. 

In  order  to  examine  whether  the  recorded  Apeel  2 ground  motion  is  atypical  for  the  recorded 
soil  profile,  a ground  motion  was  calculated  for  the  soil  profile  shown  in  figure  3.2,  assuming 
that  the  profile  shown  is  underlain  by  soft  rock,  and  using  the  Rincon  Hill  outcrop  motion. 
Since  Rincon  Hill  is  approximately  79  km  from  the  source  of  the  Loma  Prieta  Earthquake,  and 
the  source  distance  of  Apeel  2 is  47  km,  the  Rincon  Hill  acceleration  was  amplified  by  a factor 
of  2,  which  is  an  average  value  obtained  from  the  Joyner-Boore  equation  (Joyner  & Boore, 
1981)  and  from  correlations  for  rock  sites  recently  developed  for  the  Loma  Prieta  earthquake. 
Figure  5.4  shows  a comparison  between  the  response  spectra  for  the  recorded  and  calculated 
ground  motion,  which  are  quite  close.  Even  though  there  is  considerable  uncertainty  about  the 
subsurface  condition  below  the  soil  profile  shown  in  figure  3.2  and  about  the  probable 
bedrock  motion,  figure  5.4  indicates  that  it  is  not  unreasonable  to  expect  this  type  of  ground 
motion  for  a type  S3  soil  profile. 

The  response  spectra  shown  in  figure  5.3  are  for  earthquake  records  with  large  source 
distances  (47  - 80  km).  However,  earthquakes  could  originate  in  the  San  Andreas  and 
Hayward  faults  at  distances  much  closer  than  that  for  the  source  of  the  Loma  Prieta  event. 
Thus  it  is  appropriate  to  compare  the  design  spectra  for  the  San  Francisco  region  with  spectra 
generated  by,  or  likely  to  be  generated  by,  the  Santa  Cruz  outcrop  motion,  which  has  a source 
distance  of  approximately  20  km.  This  is  done  in  figures  5.5  and  5.6.  It  can  be  seen  from 
these  latter  figures,  that  for  periods  less  than  1 .5s  the  design  spectra  appear  to  be  inadequate 
for  a near-source  earthquake.  They  may  be  adequate  for  bedrock  motions,  however  soil  profile 
type  SI  also  includes  stiff  soil  sites  to  a depth  of  61  m. 

Of  particular  concern  is  the  considerable  amplification  of  the  near-source  rock  outcrop  motion 
for  periods  less  than  1 .Os,  attributable  to  most  of  the  soil  profiles  studied,  even  for  relatively 
stiff  or  shallow  soil  profiles.  This  is  not  inconsistent  with  the  observation  in  the  Loma  Prieta 
Earthquake,  that  most  of  the  structural  and  lifeline  damage  occurred  in  areas  of  deeper  soil 
deposits,  while  areas  where  bedrock  is  at  a shallow  depth  remained  relatively  undamaged 
(Lew  et  al.,  1990,  3-8).  Note  that  in  figure  5.5  the  response  spectrum  for  the  Santa  Cruz 
outcrop  motion  is  reasonably  within  the  design  spectrum  for  profile  type  SI. 

5.3  Spectral  Maps  Prepared  by  USGS 

Recently  the  United  States  Geological  Survey  (USGS)  prepared  new  earthquake  risk  maps 
which  specify  predicted  response  spectra  by  two  parameters  (Algermissen  et  al.,  1991):  The 
spectral  acceleration  ordinate  for  a 0.3s  period,  which  can  be  considered  equivalent  to  Aa  as 
used  in  NEHRP  1988;  and  the  spectral  acceleration  ordinate  at  a 1.0s  period,  which  serves 
the  same  function  as  coefficient  Av  in  NEHRP  1988,  except  that  Av  is  not  defined  for  a 
specified  period.  The  two  spectral  ordinates  define  a hypothetical  spectrum  for  a type  S2  soil 
profile.  A spectral  envelope  which  according  to  the  authors  approximates  the  actual  shape  of 
the  response  spectrum  can  be  defined  as  follows: 
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where : SA  = seismic  coefficients  for  spectral  acceleration 

n = exponent  defining  the  shape  of  the  curved  part  of  the  spectrum 

The  soil  factor  S in  this  case  is  taken  as  1 for  soil  profile  type  2 and  adjusted  accordingly  for 
other  soil  profile  types.  For  the  design  spectrum  the  authors  recommend  using  n = 1 for  T 
< Is,  and  n = 2/3  for  T > Is.  The  authors  also  imply  that  n = 1 approximates  the  actual 
predicted  spectral  shape  and  n = 2/3  is  introduced  to  preserve  the  conservatism  for  long 
period  structures  which  is  implicit  in  NEHRP  1988. 

Figures  5.7  and  5.8  compare  design  spectra  for  the  San  Francisco  region  in  accordance  with 
NEHRP  1 988  with  those  recommended  by  Algermissen  et  al.  In  figure  5.8  the  curves  for  type 
S3  and  S4  soil  profiles  are  not  plotted,  because  they  fall  outside  the  range  of  the  log-log  grid. 
Coefficients  SA  were  determined  by  scaling  the  ordinates  of  the  figure  for  San  Francisco 
presented  by  the  authors.  The  map  prepared  by  the  authors  actually  shows  higher  values  for 
SA(03)  (up  to  2 G)  over  part  of  San  Francisco,  and  higher  values  of  SA(10)  (up  to  2 G)  for  the 
Santa  Cruz  area.  It  can  be  seen  that  adoption  of  the  USGS  spectra  as  design  spectra  would 
increase  the  lateral  forces  used  in  structural  design  in  the  region.  In  figures  5.9  and  5.10,  the 
USGS  spectra  for  San  Francisco  are  compared  with  the  response  spectra  for  recorded  and 
calculated  ground  motions.  In  figure  5.10  only  near  source  ground  motions  are  included  and 
the  NEHRP  spectra  are  also  shown  for  comparison.  For  the  Santa  Cruz  rock  outcrop  motion 
the  spectra  seem  very  conservative  (it  was  noted  previously  that  present  design  spectra 
reasonably  represent  that  motion).  However,  for  periods  less  than  1.5s  the  Algermissen 
spectra  fit  the  spectra  for  deeper  soil  profiles  reasonably  well.  For  periods  greater  than  1 .5s 
the  spectra  further  increase  the  already  considerable  conservatism  inherent  in  the  present 
NEHRP  design  spectra.  For  the  San  Francisco  region,  the  change  of  the  n exponent  from  1 
to  2/3  for  T > Is  seems  conservative  in  the  extreme.  However,  this  may  not  necessarily  be 
the  case  for  other  regions  in  the  U.S. 
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Figure  5. 1 : Design  Spectral  Accelerations  Specified  by  NEHRP,  1 988,  for  the  San  Francisco 
Bay  Region 
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Figure  5.2:  Design  Spectra  Specified  by  NEHRP,  1988,  for  the  San  Francisco  Bay 
Region 


72 


Acceleration,  G's 


Period,  s 


Figure  5.3:  Comparison  of  the  NEHRP  Design  Spectra  for  Acceleration  for  the 
San  Francisco  Region  with  Response  Spectra  for  Recorded  and 
Calculated  Far  Source  Ground  Motions. 
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Figure  5.4:  Comparison  of  Response  Spectra  for  the  Recorded  and  Calculated 
Ground  Motions  at  the  Apeel  2 array. 
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Figure  5.5:  Comparison  of  the  NEHRP  Design  Spectra  for  Acceleration  for  the 
San  Francisco  Region  with  Response  Spectra  for  Recorded  and 
Calculated  Near  Source  Ground  Motions. 
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Figure  5.6:  Comparison  of  the  NEHRP  Design  Spectra  for  Acceleration  for  the 
San  Francisco  Region  with  Response  Spectra  for  Calculated  Near  Source  Ground 
Motions. 
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Figure  5.7:  Comparison  of  Design  Spectra  for  Acceleration  for  the  San  Francisco 
Region  Recommended  by  Algermissen  et  al.  1991,  with  those 
recommended  in  NEHRP,  1988. 
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Figure  5.8:  Comparison  of  Design  Spectra  for  the  San  Francisco  Region 
Recommended  by  Algermissen  et  al.  1991,  with  those  recommended 
in  NEHRP,  1988. 
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Figure  5.9:  Comparison  of  the  USGS  Design  Spectra  for  Acceleration  for  the  San 
Francisco  Region  with  Recorded  and  Calculated  Response  Spectra  for 
Far  and  Near  Source  Earthquake  Motions. 
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Figure  5.10:  Comparison  of  the  NEHRP  and  USGS  Design  Spectra  for 
Acceleration  for  the  San  Francisco  Region  with  Response  Spectra  for 
Far  Source  Earthquake  Motions. 
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6. 


SUMMARY  OF  FINDINGS 


6.1  Modeling  of  the  Ground  Motion 

A revised  version  of  the  SHAKE  program,  prepared  as  part  of  this  study,  was  used  to  calculate 
ground  surface  motions  at  the  Oakland  Outer  Harbor  Wharf  and  the  Apeel  2 USGS  array, 
where  strong  motions  caused  by  the  Loma  Prieta  earthquake  where  recorded.  In  the 
calculations  for  the  Oakland  Outer  Harbor  Wharf,  preliminary  subsurface  data  were  used,  and 
it  was  assumed,  that  the  bedrock  motion  corresponded  to  the  recorded  Yorba  Buena  rock 
outcrop  motion.  The  calculations  reasonably  predicted  the  frequency  content  of  the  recorded 
ground  motion,  but  underestimated  the  amplitude  of  the  ground  motion.  Adjustment  of  the 
strain  level  at  which  dynamic  soil  properties  are  calculated,  as  well  as  the  use  of  a range  of 
dynamic  soil  properties  did  not  materially  affect  the  response  amplitude.  Increase  in  the 
assumed  base  rock  stiffness,  for  which  no  field  data  were  available,  increased  the  amplitude 
of  the  calculated  motion.  But  only  when  the  base  rock  was  assumed  to  be  infinitely  stiff,  a 
hypothetical  condition  which  gives  an  upper  limit  for  the  base  rock  stiffness  effect  but  which 
cannot  exist,  did  the  amplitude  of  the  calculated  motion  equal  that  of  the  recorded  motion. 
Response  spectra  for  a calculated  ground  motion  at  Apeel  2,  using  the  Rincon  Hill  outcrop 
motion,  amplified  to  account  for  source  distance,  were  close  to  those  obtained  for  the 
recorded  motion.  However,  there  is  considerable  uncertainty  with  regard  to  the  subsurface 
conditions  and  the  bedrock  motion. 

The  capability  of  the  revised  SHAKE  program  to  split  the  soil  profile  into  50  layers  was  used 
to  determine  the  effect  of  reducing  the  thickness  of  soil  layers  used  in  the  analysis  to  1/4  the 
wave  length  associated  with  the  highest  frequency  component  of  the  input  motion  utilized  in 
the  calculations. 

Comparison  of  amplification  functions  calculated  for  soil  profiles  split  into  thin  layers  and 
identical  profiles  split  into  much  thicker  layers  indicated  that  the  layer  thickness  used  in  the 
analysis  had  no  noticeable  effect  on  the  results  of  the  calculations. 

6.2  Effect  of  Soil  Profile  Characteristics 

The  effect  of  various  soil  profile  types  on  the  calculated  ground  surface  motion  was  examined, 
using  input  motions  at  the  base  of  the  hypothetical  profiles  corresponding  to  the  larger 
component  of  the  horizontal  rock  outcrop  motions  recorded  at  the  Rincon  Hill,  Yorba  Buena 
and  Santa  Cruz  strong  motion  stations.  The  recorded  Rincon  Hill  and  Yorba  Buena  motions 
were  thought  to  be  characteristic  of  a far  source  earthquake,  and  the  Santa  Cruz  motion  of 
a near  source  earthquake. 

The  effect  of  the  far  source  earthquakes  on  the  various  profile  types  was  as  anticipated,  with 
strong  amplification  in  the  vicinity  of  the  characteristic  period  of  the  deposit.  Addition  of  a 24 
m layer  of  alluvium  at  the  base  of  a 30  m deposit  of  New  Bay  mud  did  not  significantly  affect 
the  predicted  ground  motion  at  the  surface  of  the  New  Bay  mud  deposit. 

All  the  soil  profiles  studied,  but  particularly  the  alluvial  deposits,  substantially  amplified  the 
near  source  outcrop  motion  for  periods  below  1 .2s.  For  periods  greater  than  1 .3s  there  was 
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no  substantial  difference  between  the  response  spectra  for  deep  New  Bay  mud  deposits 
calculated  for  the  Santa  Cruz  motion,  and  those  calculated  for  the  Rincon  Hill  motion. 

6.3  Comparison  With  Recommended  Design  Spectra 

Comparison  of  response  spectra  for  observed  and  calculated  ground  motions  with  the  NEHRP 
design  spectra  indicates  that,  for  a far  source  earthquake,  the  response  spectra  for  observed 
ground  motions  at  Oakland  Outer  Harbor  Wharf  and  Apeel  2 significantly  exceed  the 
envelopes  for  type  S2  and  S3  profiles  proposed  in  the  NEHRP  '88  recommended  design 
spectra,  as  recently  amended  in  committee. 

For  a near  source  earthquake,  it  appears  that  at  least  the  Santa  Cruz  outcrop  motion  is 
reasonably  within  the  envelope  of  the  recommended  design  spectra.  However,  the  spectra 
seem  very  unconservative  for  all  calculated  ground  motions  for  the  soil  profiles  studied.  This 
is  not  inconsistent  with  the  observation  that  in  the  Loma  Prieta  earthquake  damage  to 
buildings  and  lifelines  was  almost  exclusively  concentrated  in  areas  of  deeper  soil  deposits 
(Lew  et  al.,  1 990).  While  some  of  this  latter  damage  can  be  attributable  to  liquefaction  (not 
a topic  of  this  study),  much  of  it  occurred  in  areas  where  liquefaction  was  not  a contributing 
factor. 

All  the  response  spectra  for  calculated  and  observed  ground  motions  are  within  the  envelope 
of  the  recommended  NEHRP  '88  design  spectra  for  periods  greater  than  1.5s. 

Design  spectra  recently  recommended  in  seismic  risk  maps  prepared  by  USGS,  as  applied  to 
the  San  Francisco  Bay  region,  seem  to  reasonably  predict  the  response  spectra  for  near  source 
earthquakes  associated  with  the  soil  profiles  studied  for  periods  up  to  1.0s.  They  seem 
conservative  for  rock  outcrop  motions,  and  conservative  in  the  extreme  in  the  longer  period 
range.  It  is  not  clear  whether  these  latter  findings  can  be  extrapolated  to  other  regions  in  the 
U.S. 
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APPENDIX:  EQUATIONS  USED  IN  THE  SHAKE  PROGRAM 


Equations  used  in  the  SHAKE  program  were  derived  by  Schnabel  et  al.,  1972,  and  amended 
by  Udaka  and  Lysmer,  1973. 

For  the  horizontal  displacement: 


u=u(z,t) 


the  1-D  wave  equation  must  be  satisfied: 


d2u  n d2u 
p =G — - 

dt 2 dz 2 


+ 7] 


d3u 

dz2dt 


...  (9) 


. . . (10) 


where : u = horizontal  dynamic  displacement 

z = vertical  position  (in  direction  of  wave  propagation) 
t = time 

p = mass  density 
r ] = viscosity 

For  a simple  harmonic  displacement  with  the  circular  frequency  cj: 


u(z,t)  = U(z)-ei“[ 


aw 


where:  U = displacement  amplitude 

from  (10)  and  (11): 

(G*iuu)^f=pu2U  ...(12) 

dz 2 


The  general  solution  to  Eq(12)  is: 


U(z)  =Eelkx  +Fe  lkx 


. . . (13) 
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where : 


. . . (14) 


pcj2  _ par 
G + i cot]  G * 


k = complex  wave  number 
(j  = complex  shear  modulus 


The  critical  damping  ratio  fi  (see  figure  1 ) is  related  to  the  viscosity,  and  G*  has  been  defined 
by  Udaka  and  Lysmer  as: 


. . . (15) 


G * =G[1  -2/32  + 2//3(l  -/32)0-5] 


Within  the  frequency  range  under  consideration  G and  fi  are  assumed  to  be  independent  of 
frequency  (refer  to  figure  1).  Thus  G*  is  frequency  independent.  From  Eqs(1 1)  and  (13)  the 
wave  equation  solution  for  a harmonic  motion  of  frequency  uj  can  be  derived: 


. . . (16) 


u(z,t)  = Ee l(-kz+ut)  + Fe  z(*z 


where  E and  F are  displacement  amplitudes. 

The  first  term  in  (16)  represents  an  upward  (incident)  travelling  wave,  and  the  second  term 
a downward  (reflected)  travelling  wave.  If  the  downward  direction  of  wave  travel  is  taken  as 
positive,  then,  for  layer  "mn  of  thickness  h,  the  relative  displacement  at  the  top  and  bottom 
of  the  layer  are: 


a(0)=(£  +F)eiut 

x ' N m my 


. . . (18) 


The  shear  stress  in  a horizontal  plane  is: 
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r(z,t)  = G ’ — = ikG  ' (Ee tkx -Fe  -,tx)  e 
dz 

and  at  the  top  and  bottom  of  layer  m: 

■ • • (19) 

r(0)  =ik  G*  (E  -F  )eiot 

x 7 m m v m m/ 

. . . (20) 

r(h)=ik  G‘  (E  eik-"-+F  e e '*» 

v/  m m x m m ' 

• • • (21) 

If  stresses  and  displacements  at  layer  boundaries  are  continuous,  then: 

p j— i r~  ik  h p —ikh 

E + F =Eemm  + Fe 

m + 1 m + 1 /w  //i 

• • • (22) 

E ,-F  ,=  kmGm  (E  e,k-h~-F  e'ik-h-) 

m + 1 m + 1 * v m m ' 

+ 1 + 1 

. . . (23) 

From  Eqs.  (22)  and  (23)  recursion  formulas  for  amplitudes  E and  F can  be  derived  for 
successive  layers,  starting  with  the  top  layer: 


E =0.5  E (1+a  ) e '*”*■  + 0.5  F (1 -a  )e'ik-K 

m + 1 m'  my  m/ 

• • • (24) 

F . = 0.5E  (1  -aJeik-h-  + G.5F(\  +a)e~il'K 

. . . (25) 
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where: 


am  = frequency  independent  complex  impedance  ratio: 


otm  = 

m 


kg: 

k G * 

Km  + 1 + 1 


P G,n 
ym  m 

Pm  + 1 + 1 


0.5 


...  (26) 


On  the  free  surface  (m  = *\)  the  amplitudes  of  the  incident  and  reflected  waves  are  equal. 
Since  the  recursion  equations  (24  and  25)  are  directly  proportional  to  E and  F,  the  relative 
amplitude  of  any  layers  can  be  calculated  by  assuming  that  e1  — fl  = 1.  Starting  from  layer 
2,  em(oj)  andfm(o))  can  be  calculated  for  all  layers. 

Transfer  functions  between  the  displacements  of  any  two  layers  can  thus  be  calculated: 


_ +/„(<*) 

Un  +/» 


(27) 


Accelerations  and  strains  in  any  layer  can  be  obtained  from  displacements: 


«(z,f)  = -w  2(Ee  '<**♦»'>  =Fe 


. . . (28) 


y=ik(Ee  -Fe 


■ • - (29) 


The  elasticity  and  confinement  of  the  rock  base  is  taken  into  consideration  by  assuming  that 
the  amplitude  of  the  incident  wave  in  the  halfspace  is  independent  of  the  properties  of  the 
overlying  system.  This  assumption  is  applicable,  if  it  is  also  assumed  that  the  reflected  wave 
is  completely  absorbed  in  the  halfspace  and  therefore  does  not  contribute  to  the  incident 
wave.  If  the  halfspace  is  Layer  N,  and  N'  is  an  outcropping  rock  layer  for  which  an 
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acceleration  record  is  available,  then: 


) = = 

" V 2e„(a>) 


. . . (30) 


. . . (31) 


Equations  (30)  and  (31 ) permit  consideration  of  the  halfspace  elasticity.  If,  on  the  other  hand, 
it  is  assumed  that  the  halfspace  is  infinitely  rigid  (an  upper  limit  for  response  amplitude),  a 
rock  outcrop  motion  can  be  applied  directly  to  the  halfspace  (this  option  is  available  in  the 
program). 

The  latter  derivation  applies  to  the  case  of  steady  state  harmonic  displacements  of  circular 
frequency  oi.  In  the  program,  the  ground  motion  is  represented  by  a Fourier  transform,  and 
the  recursion  formulas,  which  are  frequency  dependent,  are  applied  to  each  term  of  the 
transform.  The  inverse  of  these  calculated  terms  represents  the  calculated  ground  motion. 

The  most  important  determination  is  the  selection  of  the  shear  modulus  and  damping  ratio. 
This  is  done  for  each  layer  for  a strain  level,  specified  as  a fraction  of  the  maximum  strain  in 
the  layer  produced  by  the  time  history.  This  requires  iterative  calculations,  which  are 
continued  for  a user  specified  number  of  iterations.  The  initial  run  is  made  for  and  a 
specified  damping  ratio.  In  subsequent  iterations  dynamic  soil  properties  are  chosen  from  sets 
of  specified  soil-specific  dynamic  soil  properties  such  as  those  shown  in  figure  1 . Thanks  to 
the  efficiency  of  the  Fast  Fourier  routine  (Cooley  and  Tukey,  1965)  and  the  fact  that  the 
curves  in  figure  1 are  very  flat  (strains  are  on  a logarithmic  scale)  and  therefore  convergence 
is  rapid  (convergence  to  within  2%  of  the  correct  value  is  usually  achieved  in  5 iterations),  the 
routine  is  efficient  and  not  very  time  consuming. 
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